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Abstract: A design for an ultra-high Q photonic crystal nanocavity 
engineered to interact with nitrogen-vacancy (NV) centers located near 
the surface of a single crystal diamond sample is presented. The structure 
is based upon a nanowire photonic crystal geometry, and consists of a 
patterned high refractive index thin film, such as gallium phosphide (GaP), 
supported by a diamond substrate. The nanocavity supports a mode with 
quality factor Q> 1.5 x 10^ and mode volume V < 0.52(A,/nGap)^> and 
promises to allow Purcell enhanced collection of spontaneous emission 
from an NV located more than 50 nm below the diamond surface. The 
nanowire photonic crystal waveguide can be used to efficiently couple light 
into and out of the cavity, or as an efficient broadband collector of NV 
phonon sideband emission. The proposed structures can be fabricated using 
existing materials and processing techniques. 
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1. Introduction 

Nitrogen-vacancy (N V) centers found in diamond are a promising system for realizing optically 
addressable solid-state spin qubits. In single-crystal diamond, negatively charged NV^ centers 
have been used in a number of experimental demonstrations relevant to quantum information 
processing, including single photon generation |[Tl|2][3l, coherent population trapping Q, and 
optical readout and manipulation of single nuclear spins ||5j |6l |2l- In order to utilize these 
properties in quantum information processing applications HHH, efficient and scalable optical 
coupling between NVs and photonic devices such as waveguides and microcavities ifTOl [TTl is 
necessary. Photonic crystal nanocavities lfT2l[T3l[T4l[T5l[T6l confine photons to sub-wavelength 
mode volumes, V, enabling Purcell enhanced coupling between an optical dipole, such as an 
NV-center or quantum dot, and the cavity mode iflTl [Tsl [191 . The coherent dipole cavity cou- 
pling rate scales as 1 / \/V, and small V is particularly beneficial in cavity-QED (Quantum 



Electrodynamic) systems where intrinsic optical loss or dipole dephasing rates are large. 

Recently, optical coupling between high-Q microcavities and NVs hosted in diamond 
nanocrystals has been reported ll20l 12111221 . and devices for studying NVs hosted in nanocrys- 
talline diamond films have been fabricated Il23ll24l . However, to-date many desirable properties 
of NVs in single-crystal diamond have not yet been observed in nanocrystalline diamond. Al- 
though spontaneous lifetime-limited optical transition linewidths have been measured in single 
NVs hosted in diamond nanocrystals, they exhibit ^10 THz inhomogeneous broadening and 
have not been measured in high yield from an ensemble of candidate nanocrystals 1251 . In bulk 
single crystal diamond, NV^ optical transitions with lifetime limited linewidths ||26|, inhomo- 
geneous broadening of less than 10 GHz, and time-averaged spectral diffusion below 100 MHz 
lID have been observed. Because of the superior properties of NVs in bulk single crystal dia- 
mond, some cavity geometries using only single-crystal diamond and air have been proposed 
ll27ll28l . However, there are difficulties in fabricating photonic devices from single-crystal dia- 
mond, since vertical optical confinement within the diamond requires either a three dimensional 
etching process, or a method for fabricating thin single-crystal diamond films |[29| . 

An alternative to fabricating devices directly from single-crystal diamond is to integrate a pat- 
terned high-index optical waveguiding layer on top of a diamond substrate, from which photons 
can couple evanescently to NVs close to the diamond surface. Recently, a hybrid GaP-diamond 
material system, consisting of an optically thin GaP film attached to a diamond substrate, was 
used to demonstrate optical coupling between NVs close to the diamond surface and ridge 
waveguides patterned in the GaP film |30|. Owing to the enhanced local density of states in the 
near field of the GaP waveguide modes, relatively efficient evanescent coupling between NVs 
and the GaP waveguide can be achieved. In order to efficiently couple to individual NVs, three 
dimensional confinement provided by an optical cavity is necessary. Here we present designs 
for a GaP-diamond photonic crystal nanowire waveguide and nanocavity which can be fab- 
ricated using available materials and processing techniques. Using numerical simulations, we 
study the sensitivity of the optical properties of these devices on structural parameters. We also 
simulate the coupling of spontaneous emission from a broadband source, such as phonon side- 
band emission from an NV-center, into the proposed waveguide design, and show that efficient 
broadband collection of NV emission is possible with these structures. 

2. Hybrid GaP-Diamond photonic crystal nanocavity 

Photonic crystal nanocavities, formed by introducing localized perturbations to planar periodic 
structures, can support ultra-high QjV resonances. Ultra-high QjY devices have been demon- 
strated at near-IR wavelengths in free standing membranes such as Si with 2 > 2.5 x 10^ lfT6l . 
In Si films supported by low-index Si02 substrates (SOI), devices with g > 1.5 x 10^ have 
been realized ISTl . Although simulations of photonic crystal cavities fabricated from diamond 
membranes («Dia ^ 2.4) predict radiation loss limited resonances with e > 10*' ESlEHlEa, 
material loss has limited Q < 600 1241 in nanocavities fabricated from nanocrystalline mem- 
branes. Photonic crystal nanocavities fabricated from single crystal diamond membranes l29l 
have not yet been demonstrated. 

Here we study photonic crystal cavities formed in a thin waveguiding layer supported by a 
diamond substrate hosting high quality NVs near the diamond surface, as indicated in Fig.[T] In 
choosing a film from which to form the waveguide layer, we are limited to materials which are 
transparent at the NV^ zero phonon transition wavelength, 'k.j^y- = 637nm, and whose refrac- 
tive index exceeds «Dia- Epitaxially grown GaP films have a nominally high refractive index 
("GaP ^ 3.3 > «Dia) and low optical absorption at 'k.-f^^/-- Photonic crystal cavities formed in 
GaP membranes have recently been demonstrated with Q ^ 1700 l33l . and in Ref. l30l, optical 
loss of ~ 72 dB/cm was measured in a GaP waveguide supported by a single crystal diamond 




Fig. 1. Schematic of tlie GaP-on-diamond phiotonic crystal cavity design, (a) Isomet- 
ric view, (b) end view. For thie optimized structure studied in this paper, [w,d,h] = 
[192, 128,640] nm, [floi^c] = [160, 141] nm (tapered over 6 periods), and the hole radius 
r = 43 nm. 



substrate. The waveguide loss measurements presented in Ref. ||30ll indicate that devices with 
absorption and surface scattering limited Q> 2x lO'* can be realized in this system. With im- 
proved processing to reduce scattering loss, it is expected that material loss limited Q exceeding 
this value should be possible. 

Realizing a high-Q photonic crystal nanocavity in this GaP-diamond system is complicated 
by the diamond substrate's moderately high refractive index (n^ia ~ 2.4). Compared to photonic 
crystal membrane or SOI devices, a diamond substrate expands the light cone into which pho- 
tons can radiate out of the nanocavity. In addition, the broken vertical mirror symmetry from the 
diamond substrate precludes the existence of a bound mode in the GaP waveguide layer for arbi- 
trarily low frequencies. As discussed below, the underlying GaP-diamond nanowire waveguide 
structure from which the nanocavity studied here is realized does not support a non-radiating 
mode in the wavelength range of interest. This is a result of the subwavelength dimensions of 
the waveguide cross-section, which provides strong modal confinement at the expense of reduc- 
ing the effective index of the waveguide modes below «Dia- However, by extending the waveg- 
uide sidewalls into the diamond substrate, radiation into the diamond substrate can be made 
arbitrarily small. As we will show below, for a realistic diamond sidewall height, the proposed 
structure supports waveguide modes whose radiation loss is smaller than the expected intrinsic 
material or scattering loss, and forms a low-loss structure for realizing a high-g nanocavity. 

2.1. GaP-diamond photonic crystal nanowire waveguide 

The underlying structure of the photonic crystal nanocavity illustrated in Fig. [T] consists of a 
GaP photonic crystal nanowire waveguide whose sidewalls have been extended into a diamond 
substrate. This waveguide supports Bloch modes which form the basis for localized resonances 
formed when the periodic symmetry of the waveguide is broken. In general, the Bloch modes 
are either leaky or guided, depending on their frequency, co, and wavenumber, k, relative to the 
air and diamond lightlines 1341 [35l . Waveguide modes with an effective index smaller than that 
of diamond, «eff = k/o3{k) < «Dia (setting c = 1), will leak into the diamond substrate. High-g 
cavity resonances are realized by engineering a perturbation which couples guided or low-loss 
Bloch modes If36ll37l . Using 3-dimensional finite difference time domain (FDTD) simulations 
ll38l to study the photonic crystal waveguide loss, we can predict an upper limit on the Q of a 
nanocavity formed from this structure. 
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Fig. 2. GaP-diamond photonic crystal waveguide band structure, and GaP, diamond sub- 
strate, and air lightlines. Shaded regions indicate the presence of a continuum of lossy 
radiating modes. In both (a) and (b) the photonic crystal waveguide modes (red and blue 
points) were calculated for a structure with h = 640 nm. The diamond lightline in (a) is 
that of a bulk diamond substrate. The filled black points in (b) are the lowest frequency Ey 
polarized modes of an infinitely tall diamond slab of width w patterned with air holes along 
the vertical axis whose spacing and radius are equal to that of the GaP waveguiding layer. 



Figure |2a) shows the first Brillouin zone of the band-structure for the lowest frequency 
GaP-diamond photonic crystal nanowire waveguide modes. Also shown are the air, diamond, 
and GaP lightlines. The calculation used waveguide cross-section dimension [w,d] = [192, 128] 
nm, hole spacing a = 160 nm, and hole radius r = 43 nm. The hole spacing was chosen to ap- 
proximately position the lowest energy waveguide mode band edge in the wavelength range of 
interest. The waveguide cross-section was not systematically optimized, and when normalized 
by the hole spacing a, is similar to that used in previous nanowire based cavity designs ifsTllSQl . 
r/a was chosen based upon initial simulations of the Q of the nanocavity design studied below 
in Sec. 12.21 and is also similar to values used in other work. A mesh with resolution a /20, per- 
fectly matched layer (PML) boundary conditions in the y and z directions (320 nm PML layer 
thickness), and Bloch boundary conditions in the x direction were used in these simulations. For 
clarity, only modes with even parity in the y direction are calculated. Odd modes lie higher in 
frequency, and as will become clear below, we are primarily interested in the lowest frequency 
waveguide mode. Note that parity in the z dimension is not conserved for these GaP-diamond 
structures, since they are not vertically symmetric. Both the lowest frequency TE-like (dom- 
inantly Ey polarized) and TM-like (dominantly Hy polarized) modes are shown in Fig. |2] We 
label these modes TE-l and TM-1, respectively. Also evident in Fig.|2tb) is a mode (black open 
points) whose dispersion lies between the TE-l and TM-1 bands. This mode corresponds to a 
leaky mode whose field is predominantly confined within the diamond substrate. Leaky modes 
with waveguide quality factor, Qwg < 100 are not shown. Qwg is related to the waveguide loss 
per unit length, a, by a = K>/Qvg, where v,^, = da/dk is the group velocity of the waveguide 
mode. 

From Fig. 12 a), it is immediately clear that for the specified waveguide dimensions, no 
waveguide modes lie below the diamond lightline. In particular, at the band edge {k = kx =K/a) 
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the effective index of the lowest frequency waveguide mode is smaller than npia; this mode is 
expected to leak into the diamond substrate. This leakage could be reduced by increasing the 
GaP waveguide dimensions (w, d), increasing the effective index of the mode at the expense of 
decreasing the peak single photon field strength, and decreasing the frequency and momentum 
mismatch between the TE-1 and higher order waveguide modes. Alternatively, we can extend 
the GaP sidewalls (including the holes) by depth h into the diamond substrate. In the limit 
that h ^ "k/iioia, the bulk diamond lightline is no longer relevant in determining whether a 
GaP-diamond waveguide mode is bound or leaky. Instead, the relevant substrate is the approx- 
imately two dimensional vertical diamond photonic crystal slab underneath the GaP photonic 
crystal waveguide. The band-structure of the lowest frequency Ey polarized mode supported by 
this vertical slab is shown in Fig. |2b). Above this band exists a continuum of radiating slab 
modes with non-zero vertical momentum; this band forms a renormalized "structured light- 
line". At k = kx, the TE-1 mode lies below this structured lightline, indicating that this mode 
is not leaky in the h S> A,/nDia limit. The existence of this mode below the structured lightline 
is not sensitive to the choice of w, since decreasing w decreases the effective index of both the 
TE-1 mode and the lowest energy radiating slab mode. Decreasing d decreases the TE-1 effec- 
tive index in a manner typical of slab and ridge waveguides ll40l . and can push the TE-1 mode 
above the structured lightline at a given frequency. 

In Fig. [3] the waveguide loss per optical cycle, l/Qwg, of the TE-1 mode is shown as a 
function of the etch depth li, at k ~ kx, for a = Oc = 141 nm, where is the "cavity" hole 
spacing of the nanocavity design presented in the following section and indicated in Fig.[T] This 
was calculated using FDTD simulations, by monitoring the power radiated into the absorbing 
y and z boundaries of a waveguide Bloch unit cell. We expect Q^g atk = kx to place an upper 
limit on the Q of a cavity mode formed predominantly from the TE- 1 mode ll37l . Figure [3] 
indicates that for h < 100 nm, the TE-1 waveguide mode is very leaky with gwg < 500, but that 
for increased etch depth h = 600 nm, the loss can be reduced and gwg > 10^. For /? > 100 nm, 
Qwg increases exponentially with h. 



TE-1 conduction band edge, a - 

TE-1 valence band edge, a - ao 

• TE-1 valence band edge, a - a^avM 
Defect mode frequency, 
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Fig. 4. (a) Lattice constant cieav as a function of position in the cavity. Each point, x,-, 
on the graph corresponds to a position midway between the center of two holes spaced 
by acavi^i)- (b) Frequency of the photonic crystal waveguide TE-1 valence band edge, 
coj^"' = (B^^''(fcx = It/a), for a set by the cavity's varying "local" hole spacing, Ocavi^i), 
shown in (a). 



2.2. Nanocavity design 

The nanocavity studied here is formed in the photonic crystal nanowire waveguide described in 
Sec. |2.1| bv locally perturbing the hole spacing, Ucav, as indicated in Figs.|4|a) and|5|a), forming 
a heterostructure BTI cavity ifTSl l42l . The hole spacing is varied slowly over 6 periods from 
Qc = 141 nm in the center of the cavity, to a constant "bulk" value of ~ 160 nm. The cavity 
is symmetric about the x and y axes. For simplicity, the hole radius remains constant. This type 
of photonic crystal nanowire waveguide cavity 1431 1441 supports localized cavity modes with 
simulated radiation loss limited Q> IQ' when applied to structures with air undercladding 
(waveguide-substrate index contrast An = «siN — «air ~ 2.0 — 1.0) B31 l46l . In SOI material 
systems (A« = nsi — "SiO? ^3.48—1 .44), nanowire based cavities with Q> Ax 10^ have been 
theoretically studied ||39ll . Ultrahigh Q/V nanocavities formed in photonic crystal nanowire 
waveguides on a high-index substrate such as diamond have not been previously investigated. 

Adjusting the local hole spacing within the photonic crystal waveguide shifts the frequency 
range of the photonic crystal waveguide stop-band. In a quasi- ID picture, the photonic crystal 
waveguide forms an "optical potential" li47J whose band edge is locally modulated by the vari- 
ation in the hole spacing or hole size Il48ll49ll46l . In the cavity design considered here, the local 
reduction in hole spacing shifts the valence and conduction band edges up in frequency so that 
the cavity supports an "acceptor" defect mode [34] formed by superpositions of the valence 
band states of the unperturbed photonic crystal nanowire waveguide 1371 . 

Following an analysis similar to that in Ref. l46l . the photonic crystal cavity optical potential 
is illustrated by Fig.|4|b), which shows the frequency CO^^"' of the TE-1 valence band edge as 
a function of the "local" hole spacing, Qcuv, at positions Xi midway between two holes in the 
cavity. The band edge frequency CO™'* was calculated at each x, by using FDTD simulations to 
determine the band structure of the nanowire photonic crystal waveguide in Sec. |2.1[ with the 
hole spacing set to acav{xi). Also shown are the frequencies of the TE-1 valence and conduction 
band edges in the "bulk" waveguide region, where a ~ «„, and of a high-Q localized cavity 
mode discussed in detail below. Note that these photonic crystal nanowire cavities have also 
been analyzed in the context of the unperturbed photonic crystal waveguide sections forming 
mirrors for photons trapped in the central cavity region l50l ISTl [39l l45l . 

FDTD simulations of this structure predict that it supports a mode at a resonance wavelength 
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Fig. 5. Dominant electric field component (Ey) of the high-2 nanocavity mode, (a) Top 
view: x — y plane bisects the GaP waveguiding layer, (b) End view: y — z plane through the 
center of the cavity, (c) Vertical profile of the field and dielectric constant. 



COc close to the NV zero phonon line at 637nm (cOc = 0.249 x 27t/fl„), with g = 1 -5 x 10^ and 
mode volume V ~ 0.52 (A,/«Gap)^, where V is defined by the peak electric field energy density: 
V — J n^\E\-^dr/ (n^jZip)^^^. The maximum field amplitude inside the diamond is Es = 0.72Eo, 
where £„ is the field maximum at the location of maximum energy density. E„ is located inside 
the high index GaP layer for the device and mode studied here. The field maximum decays 
to O.lOiio at a depth z = 155 nm below the diamond surface. Simulations indicate that similar 
designs with the number of grading periods reduced to 4 support more tightly confined modes 
(V ^ 0.3 (A-/«Gap)^) with a maximum Q ^ 6 x 10"*. The mesh used in these simulations had 
uniform resolution fl„/20, and extended = 8 periods («„) in ±i beyond the end of the graded 
region, SOf, in ±£ above and below the diamond surface, and 3flo along the ±y directions. A 
PML layer thickness of A-Og was used in the highest-g simulations. After exciting the structure 
with a narrow band source, the steady state ratio of stored and radiated power was used to 
calculate the modal Q. The radiated power along each axis was monitored using reference 
frames positioned ao/2 away from each PML boundary. Mirror symmetry was enforced about 
the X (odd parity) and y (even parity) axes. Higher resolution (ao/25) simulations of the h = 640 
nm structures were also conducted, and it was verified that the structure supports a mode with 
e> 1.5x 10^\ 

The role of the diamond substrate in limiting Q is illustrated by simulating its dependence 
on the depth of the etched diamond ridge. Figure|6|a) shows 2 as a function of li, and indicates 
the relative contributions to Q due to radiation into the bottom diamond substrate {Q\,ot), top 
air cladding (Qtop), end of the cavity (Qend) and side of the cavity (Qside)- The total cavity Q is 
given by g ' = Gtop + Gbot + Qend + 2side- When h is reduced to 160 nm, increased radiation 
loss into the substrate degrades Q < 500. For h > 640 nm, radiation loss through the end of the 
cavity into the photonic crystal waveguide begins to play a non-negligible role in limiting Q. 

Figure |6tb) shows the dependence of Q on the number of waveguide periods between 
the cavity grading edge and the PML absorbing layer used in the simulation. For > 8, Qend 
is approximately constant, indicating that it is limited by coupling between the cavity mode 
and either photonic crystal waveguide modes which do not exhibit a bandgap at co^ (e g-, the 
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Fig. 6. Dependence of nanocavity mode Q on geometric parameters. The total Q (blue 
points) and the contribution to Q from radiation into specific directions (other colored 
points) are given. Dependence of Q on (a) etch depth h, (b) number of "mirror" periods 
N in the x direction between the edge of the graded cavity region and the PML absorbing 
simulation boundary, (c) cavity minimum (center) hole spacing Oc- For A' < in (b), the 
PML boundary of the simulation domain overlaps |A'| periods of the graded cavity region. 



TM-1 mode) or radiation modes ISOllSTI . For N < 5, the cavity loss is dominated by "mirror" 
leakage. In a practical device, represents the number of "mirror" periods between the edge 
of the cavity grading and an integrated waveguide, and would be set to a value that ensures that 
radiation loss from of the cavity is predominantly into the waveguide. 

The sensitivity of Q to the central cavity hole spacing, Oc, is shown in Fig. ^c). When ac 
is varied between 132— 150 nm, Q remains > 5 x 10^. In these simulations, the bulk hole 
spacing is maintained fixed at Oo = 160 nm, and the grading of Qc remains parabolic. Q and 
2end are maximized when = 141 nm, while Qbot trends upward as Uc — > ag. In this limit, 
the depth of the parabolic cavity defect becomes increasingly small, minimizing coupling to 
lossy waveguide modes close to the diamond lightline. As a result, Qbot increases. Qend is 
maximized when a sufficiently small is chosen to position (0^ well within the TE-1 bandgap, 
as illustrated in Fig.|4|b). However, pushing co^ too far above the TE-1 valence band edge can 
result in efficient phase matched coupling between the cavity mode and the lossy TM-1 valence 
band edge modes (see Fig. |2|, reducing Qend- These competing dependencies determine the 



optimum Gc- Note that V varies between 0.42 — 0.70 (^/noap)^ over ttie simulated range of Uc- 



2.3. Nanocavity cavity-QED parameters for coupling to NV centers 

From the FDTD calculations of V and Q presented above, we can predict the coherent coupling 
rate, g, for a single NV^ center located near the cavity and coupled to a single photon stored 
in the nanocavity mode. For experiments in cavity QED we must be careful to distinguish be- 
tween zero-phonon optical transitions (these become spectrally narrow at low temperature) and 
phonon-assisted transitions. The total spontaneous emission rate of an NV^ center is measured 
to be Ytot = 27t X 13 MHz; however the rate into the ZPL alone, Yzpl, is only 3% of Ytot 1521 . 
The coherent coupling rate between a single photon and the NV^ ZPL is given by 
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where V = V / (k / ncap)^ , |£'(rNv)| is the magnitude of the electric field at the NV location, and 
Eo is the field strength at the electric field energy density maximum, which is located within the 
GaP for the device studied here. This expression assumes that the spatial orientation of the NV 
transition dipole moment and electric field polarization are parallel, which is possible for one 
of four allowed NV orientations in (111) diamond. For the cavity mode volume given in Sec. 
12.21 Eq. ([T]) gives gs/ln = 2.25 GHz for an NV optimally located at the diamond surface. 

Assuming that the cavity Q is limited by radiation loss as calculated in Sec. 12.21 we see 
that [gj,K,ytot] /27t = [2.25,0.16,0.013] GHz, where K = Co/2g, indicating that it is possible 
to reach the strong coupling regime |[53l using this nanocavity design. Also of interest is the 
Purcell enhanced spontaneous emission rate of the NV^ ZPL into the cavity mode. Assuming 
the cavity and NV^ ZPL are aligned spectrally, the Purcell enhancement factor to the ZPL 
spontaneous emission rate is given by. 
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For an NV optimally positioned at the diamond surface, Eq. (|2]i gives = 4.9 x 10^. 

In a realistic experiment, it may not be possible to couple to NVs with desirable optical and 
spin coherence properties arbitrarily close to the diamond surface |54j|. In addition, material 
loss and fabrication imperfections will likely limit the nanocavity Q below the highest values 
presented in Sec. 12.21 If we assume that the NV is optimally positioned 50 nm below the dia- 
mond surface, and that the cavity Q is limited to 2 x lO'* (e.g., due to material absorption and 
surface scattering ||30l ), we find F > 16, indicating that ~ 94% of the total NV^ spontaneous 
emission will radiate into the nanocavity mode. The 6% of uncoupled NV emission will radiate 
into phonon sidebands not resonant with the cavity mode. 

2.4. Waveguide collection ofNV phonon sideband emission 

In addition to forming the basic structure for the the nanocavity presented in Sec. 12.21 the 
waveguide studied in Sec. 12.1 l ean function as an efficient broadband collector of radiation from 
NV phonon sidebands. Measurements of broadband phonon-sideband emission are used in ex- 
periments to gather information regarding the relative populations of the NV spin triplet ground 
states im, as well as the nuclear spin state of neighboring impurities ||5]|6]|7]. From coupled- 
mode theory POl , the normalized power spectral density, 15(03)1', coupled into a waveguide 
mode from a NV dipole located at position Tnv in the photonic crystal waveguide near field is 
approximately given by 
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Fig. 7. (a) Geometry of the simulated dipole and the asymmetric photonic crystal waveguide 
hole pattern. Left of the dipole, the hole spacing is set to a value (a,,) larger than the hole 
spacing (oc) to the right of the dipole. For frequencies close to the valence band edge coj^ ' 
of the right waveguide region, there are no propagating waveguide modes to the left of the 
dipole. In the transition region, the hole spacing and hole radius are graded to 0.75 of their 
nominal values in order to reduce back reflections 151| . More precisely, the hole spacings 
and hole radii in the transition region are scaled by [0.95,0.86,0.8,0.75] (moving in the 
+x direction) from the values in the waveguide region, (b) Normalized radiated spectra of 
the dipole into waveguide and non-waveguide radiation channels, (c) Normalized radiated 
spectra of the dipole into the principal directions of the simulation domain. 



where «^ (co) is the frequency dependent group index of the waveguide mode, £ (r ) is the waveg- 
uide mode field amplitude, and A is the waveguide mode area defined by 

«"(«(r)2|£(r)|2) 

\ / max 

where volume u is a unit cell of the waveguide. We have assumed that the dipole and waveguide 
field at Tnv are parallel. |5(co)p is normalized by the "bulk" power spectral density of the 
NV dipole embedded far below an unpatterned diamond surface. Equation (O indicates that 
coupling into a given waveguide mode can be increased by maximizing iig/A. Photonic crystal 
waveguides are ideally suited in this regard 155)1 : they can support modes with sub-wavelength 
modal area and large rig. 

To gain a quantitative estimate of the NV- waveguide coupling, and to study enhancement or 
suppression of N V coupling into other radiation modes of the waveguide structure, it is useful to 
conduct FDTD simulations of the emission spectrum of a broadband dipole source positioned 
at Tnv. Here we consider a broadband dipole, polarized along y and located znv = 11 nm below 
the diamond surface and midway between two holes of the photonic crystal waveguide studied 



in Sec. 12. II In order to channel the dipole emission into a forward propagating waveguide mode, 
the hole spacing of the waveguide structure used here is asymmetric about the dipole position, 
as shown in Fig. 13^)- For «„ > cic, the valence band edge of this structure is lower on the 
left side of the dipole than on the right side; the band edge of the left ("mirror") region has a 
lower frequency than the band edge of the right ("waveguide") region. For efficient collection 
of NV emission, the structure should be designed such that the majority of the NV phonon 
sideband emission is lower in frequency than the "waveguide" TE-1 valence band edge, and 
higher in frequency than the "mirror" TE-1 valence band edge. Note that reflections from the 
waveguide termination can dramatically alter the coupling between a dipole and the waveguide, 
as studied in ||55| . Here we have chosen a termination designed to suppress these reflections 
in the simulations. In FDTD simulations of periodic structures supporting propagating modes, 
PML terminations of the propagating axis can create significant reflections. A waveguide which 
is invariant along the propagating dimension, such as the termination region of the structure 
considered here, is better suited for a PML boundary. 

Figure |2tb) shows the FDTD calculated spectrum for the system described above. Power 
spectral density of radiation into the waveguide, and into all other channels (e.g., the substrate) 
are shown. In the simulations studied here, Uc = 141 nm and «„ = 160 nm. The hole radius 
and waveguide cross-section are the same as in Sec. 12. II and h = Aa,,. The spectra are normal- 
ized by the spectrum of the identical dipole source when it is positioned at depth znv below 
an unpatterned diamond surface. From Fig.|7tb), we see that for frequencies below the valence 
band edge of the waveguide region (co < coj^"' ~ 0.25 x 2%/ao), the emission into the waveg- 
uide exceeds the total "bulk" emission in absence of the waveguide. The radiated waveguide 
power was calculated by monitoring the power flux through an area overlapping the waveguide 
cross-section, with dimensions 2w x 2d, located near the simulation PML boundary. Note that 
emission into non-waveguide modes is also enhanced. As shown in Fig.|7jc), radiation into the 
substrate (— £ direction) and into the waveguide (+x direction) are the dominant radiation chan- 
nels. Also, note that relatively little power is coupled into the backward propagating waveguide 
mode (— i direction) in the displayed frequency range. For Ug = 160 nm, the bandwidth of this 
efficient waveguide coupling is approximately 637 — 700 nm, allowing efficient collection of a 
large portion of the NV^ zero phonon line and phonon sideband emission. This may be partic- 
ularly useful in room-temperature NV experiments which do not rely upon readout of a sharp 
optical ZPL, but require efficient collection of sideband emission, for electron spin readout, for 
example. 

2.5. Fabrication considerations 

Fabrication of the devices studied in this paper requires three key processing abilities: (i) pat- 
terning GaP or another high index film with the design presented above, (ii) attaching high index 
films to a diamond substrate, and (iii) transferring the thin film pattern into the diamond sub- 
strate. Proofs of principle of these individual steps have been demonstrated. Results presented 
in Refs. [33. 24J show that it is possible to lithographically define and etch sufficiently small 
holes using current electron beamwriting and plasma etching tools. Large area patterned GaP 
films were transferred and adhered 1561 onto a diamond sample in Ref. ||30|, as have smaller 
micron-scale structures ifSTl . Transferring the pattern into the diamond can be accomplished us- 
ing oxygen based plasma etching, e.g., in Ref. 1541 a highly selective plasma etch was used to 
extend a SiN pattern into a diamond substrate. The degree to which the precise design studied in 
Sec. l2.2l can be replicated will affect the maximum achievable Q, and achieving the high aspect 
ratio diamond sidewalls of the h = 640 nm structure may be challenging. Theoretical studies 
analyzing the importance of the sidewall slope may be necessary in future work. Finally, note 
that although precise (2 nm) tuning of ac was considered in Sec. 12.21 in order to maximize Q, 



the cavity supports modes with Q> 5 x 10^ over the range 132 nm < Oc < 150 nm. 

The GaP-based cavity proposed here is designed specifically for coupling to negatively 
charged NV centers positioned within 50 nm of a diamond surface. It is well established that a 
dense layer of NV centers close to a diamond surface can be created in high-nitrogen diamond 
by ion implantation and annealing Il58l[59ll60ll54l . Such a dense layer is well-suited for initial 
testing of the optical characteristics of GaP-based microcavities. However, for the proposed de- 
vices to be useful for quantum information applications, single NV centers with good spectral 
characteristics must be fabricated close to a surface in high-purity diamond. Initial tests suggest 
that in high-purity diamond, a nearby surface can have a large effect on NV properties 1541 that 
include charge stability and optical linewidth. A thorough study of how to optimize the proper- 
ties of single NV centers close to a surface will be required for any device based on NV centers 
coupled to an optical microcavity, including the GaP-based design presented here. 

3. Conclusion 

In this paper we have analyzed a photonic crystal geometry which supports low-loss, sub- 
wavelength nanocavities and waveguides suitable for coupling to diamond NV-centers in sin- 
gle crystal diamond. The structures presented here can be fabricated using existing materials 
and processing techniques, and promise to allow efficient collection of photons emitted from 
NV-centers, a crucial requirement for proposed applications Q . Efficient collection of N V ZPL 
emission is necessary for single photon and quantum repeater applications of N Vs, and efficient 
collection of NV sidebands is required for high bandwidth readout of the NV electron spin. The 
photonic crystal nanocavity design presented above supports modes with the necessary Q and 
V to reach the strong coupling or large Purcell factor regimes with the ZPL of an NV-center. 
The photonic crystal waveguide analyzed here allows broadband collection of NV emission, 
and is promising for efficient collection of NV sidebands. In future work, it is expected that in- 
tegration of multiple cavities, and realization of nanophotonic circuits for quantum information 
processing using NV-centers, will be possible using the device designs presented here. 



This figure "pc_cavity_iso_end.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/0904.0500v2 



